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Abstract 

Noncommutative multi-solitons are investigated in Euclidean two-dimensional U(n) and Grass- 
mannian sigma models, using the auxiliary Fock-space formalism. Their construction and mod- 
uli spaces are reviewed in some detail, unifying abelian and nonabelian configurations. The 
analysis of linear perturbations around these backgrounds reveals an unstable mode for the 
U(n) models but shows stability for the Grassmannian case. For multi-solitons which are diag- 
onal in the Fock-space basis we explicitly evaluate the spectrum of the Hessian and identify all 
zero modes. It is very suggestive but remains to be proven that our results qualitatively extend 
to the entire multi-soliton moduli space. 
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1 Introduction 



Multi-solitons in noncommutative Euclidean two-dimensional sigma models are of interest as static 
DO-branes inside D2-branes with a constant B-field background [1] but also per se as nonperturba- 
tive classical field configurations. Adding a temporal dimension, they represent static solutions not 
only of a (WZW-extended) sigma model but also of the Yang-Mills-Higgs BPS equations on non- 
commutative R 2+1 . In fact, the full BPS sector of the Yang-Mills-Hi ggs system is, in a particular 
gauge, given by the (time-dependent) solutions to the sigma-model equations of motion. 

In the commutative case, the classical solutions of Euclidean two-dimensional sigma models 
have been investigated intensively by physicists as well as by mathematicians (for a review see, 
e.g. [2]). Prominent target spaces are U(n) group manifolds or their Grassmannian cosets Gr(n, r) = 
U(r)xU(i-r) for 1 < r < n (which are geodesic submanifolds of U(n)). 1 Any classical solution <I> 
for the U(n) sigma model can be constructed iteratively, with at most n— 1 so-called unitons as 
building blocks [3, 4]. The subset of hermitian solutions, & = <I>, coincides with the space of 
Grassmannian solutions, <E> 6 Gr(n, r) for some rank r < n, of which again a subset is distinguished 
by a BPS property. 2 These BPS configurations are precisely the one-uniton solutions and can be 
interpreted as (static) multi-solitons. 

Let us now perturb such multi-solitons within the configuration space of the two-dimensional 
(i.e. static) sigma-model, either within their Grassmannian or, more widely, within the whole group 
manifold. A linear stability analysis then admits a two-fold interpretation. First, it is relevant for 
the semiclassical evaluation of the Euclidean path integral, revealing potential quantum instabilities 
of the two-dimensional model. Second, it yields the (infinitesimal) time evolution of fluctuations 
around the static multi-soliton in the time-extended three-dimensional theory, indicating classical 
instabilities if they are present. More concretely, any static perturbation of a classical configuration 
can be taken as (part of the) Cauchy data for a classical time evolution, and any negative eigenvalue 
of the quadratic fluctuation operator will give rise to an exponential runaway behavior, at least 
within the linear response regime. Furthermore, fluctuation zero modes are expected to belong to 
moduli perturbations of the classical configuration under consideration. The current knowledge on 
the effect of quantum fluctuations is summarized in [2]. 

The Moyal deformation of Euclidean two-dimensional sigma models has been described in [5], 
including the construction of their BPS solutions. The operator formulation of the noncommutative 
U(n) theory turns it into a (zero-dimensional) matrix model with U(C n <g>W) as its target space, 
where 7i denotes a Heisenberg algebra representation module. Each BPS configuration again 
belongs to some Grassmannian, whose rank r can be finite or infinite. The latter case represents a 
smooth deformation of the known commutative multi-solitons, while the former situation realizes 
a noncommutative novelty, namely abelian multi-solitons (they even occur for the U(l) model). 
Although these solutions are available in explicit form, very little is known about their stability. 3 
Our work sheds some light on this issue, in particular for the interesting abelian case. 

The paper is organized as follows. In the next section we review the noncommutative two- 
dimensional U(n) sigma-model and present nested classes of finite-energy solutions to its equation 
of motion. We focus on BPS configurations (multi-solitons) and their moduli spaces, unifying the 
previously different descriptions of abelian and nonabelian multi-solitons. Section 3 is devoted to 
a study of the fluctuations around generic multi-solitons and detects a universal unstable mode. 

^ost familiar are the Gr(n, 1) = CP™ -1 models. 

2 not to be confused with the BPS condition for the Yang-Mills-Higgs system 

3 A few recent works [6, 7, 8, 9, 10, 11] address perturbations of noncommutative solitons but not for sigma models. 
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Section 4 analyzes the spectrum of the Hessian specifically for U(l) backgrounds and settles the issue 
for diagonal BPS solutions, analytically as well as numerically. Section 5 addresses the fluctuation 
problem in the nonabelian case for the example of the noncommutative U(2) model. We close with 
a summary and a list of open questions. 

2 Noncommutative 2d sigma model and its solutions 

Before addressing fluctuations, it is necessary to present the noncommutative two-dimensional 
sigma model and its multi-soliton solutions which we will set out to perturb later on. These 
solutions have been discussed as static solutions of a 2+1 dimensional noncommutative sigma 
model [5, 12, 13, 14] which results from Moyal deforming the WZW-modified integrable sigma 
model [15, 16, 17]. Here we review the results pertaining to the static situation. 

2.1 Noncommutativity 

A Moyal deformation of Euclidean R 2 with coordinates (x, y) is achieved by replacing the ordinary 
pointwise product of smooth functions on it with the noncommutative but associative Moyal star 
product. The latter is characterized by a constant positive real parameter 9 which prominently 
appears in the star commutation relation between the coordinates, 

x-ky-y-kx = [z,y]* = i0. (2.1) 

For a concise treatment of the Moyal star product see [18, 19, 20]. It is convenient to work with 
the complex combinations 

z = x + iy and z = x — iy => [2,2]* = 28 (2-2) 

and to scale them to 

a = -fy and «t = -L [ a , a t], = 1. (2.3) 

A different realization of this Heisenberg algebra promotes the coordinates (and thus all their 
functions) to noncommuting operators acting on an auxiliary Fock space Ti but keeps the ordinary 
operator product. The Fock space is a Hilbert space with orthonormal basis states 

\m) = -4= (a f ) m |0) for m G N and a\0) = , 

, , . , ( 2 - 4 ) 

a\m) = \pm\m— 1) , cv \m) = \/m+l |m+l) , N \m) :=a'a\m) = m\m) , 

therewith characterizing a and a) as standard annihilation and creation operators. The star- 
product and operator formulations are tightly connected through the Moyal- Weyl map: Coordinate 
derivatives correspond to commutators with coordinate operators, 

V29d z -ad(a t ) , V26d 2 ad(a) , (2.5) 

and the integral over the noncommutative plane reads 

fd 2 xU{x) = 2vr0Tr w / op , (2.6) 

where the function /* corresponds to the operator / op via the Moyal- Weyl map and the trace is 
over the Fock space Ti.. We shall work with the operator formalism but refrain from introducing 
special notation indicating operators, so all objects are operator-valued if not said otherwise. 
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2.2 Two-dimensional sigma model 

The fields $ of the noncommutative two-dimensional U(n) sigma model are unitary nxn matrices 
with operator- valued entries, i.e. <1> € U(C™ <g) H) = \J(H® n ), and thus subject to the constraint 

$$t = i n ® i w = $t$ . (2.7) 

The Euclidean action of the model coincides with the energy functional of its 2+1 dimensional 
extension evaluated on static configurations [5], 4 

E[$] = 2TT0Tv(d z ^ ds^ + dz^ d z <S>) 

= 7rTr([a,$]t[a,$] + [a,$t][ a ,$t]t) 

= 2vrTr([a,$] t [a, $]) = 2vr|[a,$]| 2 
= vrT^A^^ + ^A^) , 

where the trace is taken over the Fock space 7i as well as over the U(n) group space. Here, we have 
introduced the hermitian Laplace operator —A which is defined via 

AO := [a,[a\0]] = [a\[a,0]] . (2.9) 

Its kernel is spanned by functions only of a or only of cl< . Varying E[&\ under the above constraint 
one finds the equation of motion, 5 

= [a,$V>$]] + [aWta,*]] = ^ A$ - A$ f $ . (2.10) 

Since the annihilation and creation operators above should be read as l„0a and l n (g) a) , 
respectively, this model enjoys a global U(l)xSU(n)xSU(n) invariance under 

$ — ► (V ® l n )$(W® l n ) (2.11) 

with VV^ = t n = WW^ , which generates a moduli space to any nontrivial solution of (2.10). 
Another obvious symmetry is induced by the ISO (2) Euclidean group transformations of the non- 
commutative plane, as generated by the adjoint action of a, a' and N. Specifically, a global 
translation of the noncommutative plane, 

(z,z) — ► (z + C, z + Q = VW(a + a, o f + a) , (2.12) 

induces on (operator-valued) scalar functions / the unitary operation 

/ ► e _Ca *- f8 */ = e aad(ffl,) -" ad(fl) / = e ^-aa fe - a a^+aa = . D( a ) t . (2.13) 

Its action on states leads to coherent states, e.g. 

£>(a)|0) = e aat "" a |0) = e"^ a e Qflt |0) =: |a) . (2.14) 
A global rotation of the noncommutative plane, 

(z,z) — ► (e i,? z, e" 1 ^) = >/20(e*V e _i V) , (2.15) 

4 |^4| 2 = Tr(^4^) ig the squared Hilbert-Schmidt norm of U(n)-valued operators on 7i. We restrict ourselves to 
finite-energy configurations, i.e. we demand that [a, $] exists and is Hilbert-Schmidt. Furthermore we only consider 
solutions for which A3> is traceclass and $ is a bounded operator in order for the below expressions to be well defined. 

5 assuming that the appearing total derivative terms vanish 
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induces the unitary transformation 

/ _^ e Wad(ota)^ = e Wotoy: e -Wota = . / )t ( 2 .16) 

for G R/27rZ, because [AT", /] = zd z f — d z fz. Applying R($) to a coherent state we obtain 

R(0)\a) = \e^a) . (2.17) 

Since the adjoint actions of both D(a) and R{$) commute with that of A, the energy functional 
is invariant under them, and all translates and rotations of a solution to (2.10) also qualify as 
solutions, with equal energy. However, other unitary transformations will in general change the 
value of E. 

There is a wealth of finite-energy configurations <!> which fulfil the equation of motion (2.10). 
A subclass of those is distinguished by possessing a smooth commutative limit, in which $ merges 
with a commutative solution. The latter have been classified by [3, 4]. We call these "nonabelian" 
because they cannot appear in the U(l) case, where the commutative model is a free field theory 
and does not admit finite-energy solutions. Yet, there can (and do) exist non-trivial abelian finite- 
energy solutions at finite 9, whose 6 — > limit is necessarily singular and produces a discontinuous 
configuration. We term such solutions "abelian" even in case they are imbedded in a nonabelian 
group. 

Of particular interest are configurations diagonal in the oscillator basis (2.4) as well as in C n , 
namely $ = ®" =1 diag({e ia i }^ ) w ith of € It. A short computation reveals that such a configu- 
ration obeys the equation of motion (2.10) if and only if the phases satisfy e loli = ±e ia * where the 
sign depends on I. Its energy is finite if the number of positive signs or the number of negative 
signs in each block labelled by i is finite. 

2.3 Grassmannian configurations 

It is a daunting task to classify all solutions to the full equation of motion (2.10), and we shall focus 
on the subset of hermitian ones. Any (not necessarily classical) hermitian configuration obeys 

$t = $ =^ $ 2 = t n ® t n =: 1 , (2.18) 

and is conveniently parametrized by a hermitian projector P=P^ via 

$ =: 1-2P = P ± - P = e inP with P 2 = P , (2.19) 

where P 1 - = 1— P denotes the complementary projector. This allows one to define a topological 
charge like in the commutative case as [21] 

Q[$] = \OTr(<S>d z <&dz<S>-<S>dz<S>d z <&) 

= Tr(P[a?,P][a,P]-P[a,P][a?,P]) 

= \P[a,P]\ 2 -\[a,P]P\ 2 (2.20) 

= Tr(P - PaPalP + PatPaP) 

= Tx(Pa{l-P)a ] P - Pa\t-P)aP) , 
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which may be compared with the energy 



= Tr([<J,P] [P,a]) 

= \P[a,P]\ 2 + \[a,P]P\ 2 (2.21) 

= Tr(P(aa t +o t o)P - PaPa)P - Pa)PaP) 

= Ti{Pa{l-P)a)P + Pa){l-P)aP) . 

Each hermitian projector P gives rise to the complementary projector 1 — P, with the properties 6 

E[l-P] = E[P] and Q[l-P] = -Q[P] . (2.22) 
Comparing (2.20) and (2.21) we get the relations 

±E[P] = Q[P]+2\[a,P]P\ 2 = -Q[P]+2\P[a,P]\ 2 , (2.23) 
which yield the BPS bound 

E[P] > 8tt\Q[P}\ (2.24) 

for any hermitian configuration. 

The set of all projectors unitarily equivalent to P is called the Grassmannian, and it is given 
by the coset 

Gr(P) = — > ; . (2.25) 

v ; U(imP) x U(kerP) v ; 

Thus, each given P (and thus hermitian <5) belongs to a certain Grassmannian, and the space of 
all hermitian <E> decomposes into a disjoint union of Grassmannians. The restriction to hermitian 
reduces the unitary to a Grassmannian sigma model, whose configuration space is parametrized by 
projectors P. In the commutative case the topological charge Q is an element of 7T2(Gr(P)) = Z, 
where S 2 is the compactified plane. Hence, it is an invariant of the Grassmannian if one excludes 
from (2.25) any singular unitary transformations with nontrivial winding at infinity. It is less 
obvious how to properly extend this consideration to the infinite-dimensional cases encountered 
here [22]. We therefore take a pragmatic viewpoint and demand Q[P] to be constant throughout 
the Grassmannian. This may downsize the above coset Gr(P) by restricting the set of admissible 
unitaries U. 

Let us make this more explicit by computing Q[UPW]. To this end, we define 
to := rf[a,U] and J" = t/t [a f , U] , with [a + u , a) + lo ] ] = 1, (2.26) 
as elements of the Lie algebra of \J(H® n ). 7 A short calculation yields 

Q[UPU*] = Tr(P[a f W,P] [a+u,P] -P[a+u,P] [<J+u*,P]) (2.27) 
= Q[P] + Tr P ( [a\ P]lo + u[a\ P] - [a, P]J-J[a,P] - uj\\-P)uj + u(t-P)J)P , 



By a slight abuse of notation, we denote the energy and topological charge as functionals of P again with the 
symbols E and Q, respectively. 

7 In fact, the finite-energy condition (see footnote before (2.8)) enforces [a,UPW] to be Hilbert-Schmidt which 
implies that [u , P] is Hilbert-Schmidt as well. Unfortunately, this does not suffice to guarantee the constancy of Q 
in Gr(P). On the other hand, u) itself need not even be bounded, as in the example of (2.16) where u> = (e 11 ' — l)o. 
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which constrains oj in terms of P. In case of TrP < oo, this indeed reduces to 

Q[UPrf] = Q[P] + Tr P ([u, at] + [a, J] + J]) = Q[P] . 



(2.28) 



Alternatively, we may calculate the infinitesimal variation of Q[P] under P — > P+5P. Remem- 
bering that P and 5P are bounded and that their commutators with a or a) are Hilbert-Schmidt, 

we get 

6Q[P] = Tr([<J, [P,8P][a,P]]-[a,[P,6P][<J,P]]) + 3 Tr 6P([a\ P] [a, P] - [a, P][a) , P]) (2.29) 
with the first trace being a "boundary term" . Variations inside the Grassmannian are given by 

5P = [A Q ,P] with Al = -A Q and PA Q P = = (1-P)A (1-P) (2.30) 
(see Section 3.2 below), and thus 

SQ[P] = Tr([a,A G [at,P]] - [at , A Q [a, P}]) + 3 Tr [P , A G [a, P] [a\ P] - K [a\ P] [a, P}] . (2.31) 

Hence, if A Q is bounded, the second term vanishes and SQ[P] reduces to the boundary term. 8 
Formally, another invariant is the rank of the projector, 

r := rank(P) = dim(imP) , (2.32) 

which may differ from Q when being infinite. Let us consider the class of projectors which can be 
decomposed as 



P = U (P + P') C/T with P = P <g> t n 



and 



PP' = = P'P 



(2.33) 



where U is an admissible unitary transformation, P denotes a constant projector, and the rank r' 
of P' is finite. In this case the difference of r and Q is determined by the invariant U(n) trace [23] 



R[P] = trP € {0,l,...,n} , 



(2.34) 



and such projectors are characterized by the pair (R,Q). Formally, the total rank then becomes 
r = P-oo + Q. Since a and a) commute with P, the topological charge depends only on the 
finite-rank piece, 

Q[P] = Q[P+P'} = Q[P'] = TrP' = r' . (2.35) 

For P=0 one obtains the abelian solutions because they fit into \J(H). For P>0 we have nonabelian 
solutions because more than one copy of H is needed to accomodate them. 

Any projector in C n <S>7i = 7i® n can be parametrized as 

P = \T) (r|r) _1 (T| (2.36) 



where 



\T) = (JT 1 ) \T 2 ) 



1) 



(\Tl) \T?) 

n) it|) 

\\Tl) |T n 2 ) 



\T[)\ 

m 

\T r n )) 



\T 2 ) 

\\Tn)J 



(2.37) 



3 Yet again, this condition is too strong a demand as the examples of (2.13) and (2.16) indicate. 
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denotes an nxr array of kets in H, with r possibly being infinite. Thus, 

(T\T) = (V|T m >) = (jXiPfPT)) 



(2.38) 



stands for an invertible rxr matrix, and r is the rank of P. The column vectors |T*) span the 
image imP of the projector in 7i® n . There is some ambiguity in the definition (2.36) of |T) since 



|T) 



\T) T for r G GL(r) 



(2.39) 



amounts to a change of basis in imP and does not change the projector P. This freedom may 
be used to normalize (T\T) = t r , which is still compatible with T G U(r). On the other hand, a 
unitary transformation 

|T) -► U\T) for f/£U(H ffi ") (2.40) 
yields a unitarily equivalent projector UPW. Altogether, we have the bijection 

p = UPU ] |f) = C/S|T) = C/|T)r with S G GL(imP) . (2.41) 

Here, the trivial action of U G U(imP) xU(kerP) on P can be subsumed in the S action on \T). 



Any diagonal projector can be cast into the form 

P d = diag . . . , l w , Pq,0 h ,...,0 h ) with 



|m)(m| 



(2.42) 



R times n—R—l times 

/en 



m=0 



Formally, Pj has rank r in 7^® ra . A corresponding rxn array of kets via (2.36) would be 

• • On \ 



(\H) 
\H) 



\ 



o Q 
|t ) 

o Q y 



with 





= (|0> |1) |2> 


|3) . 


••) 





= (000 


. 


••) 


\Tq) 


= (|o> |i) 


IQ- 


-i)) 




= ( .. 


. 


) 



(2.43) 



Q times 



In the abelian case, n=l and P=0, this reduces to \T&) = \Tq), and any projector is unitarily 
equivalent to Pq. Due to (2.22), analogous results are valid in the complementary case P — > 1— P. 

2.4 BPS solutions 

Let us focus on classical solutions within a Grassmannian. For <E> = 1— 2P, the equation of mo- 
tion (2.10) reduces to 

= [AP,P] = [a f , (l-P)aP] + [a, Pa f (l-P)] = [a , (l-P)a f P] + [a f , Pa(l-P)] . (2.44) 

Still, we do not know how to characterize its full solution space. However, (2.44) is identically 
satisfied by projectors subject to [24, 9] 



either the BPS equation 
or the anti-BPS equation 



= [a,P]P = (t-P)aP 
= [a\P}P = (l-P)a t P 



(2.45) 
(2.46) 
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which are only "first-order". Solutions to (2.45) are called solitons while those to (2.46) are named 
anti-solitons. Hermitian conjugation shows that the latter are obtained from the former by ex- 
changing P <-> 1— P, and so we can ignore the anti-BPS solutions for most of the paper. Equation 
(2.45) means that 

a maps imP imP (2.47) 

and, hence, characterizes subspaces of H® n which are stable under the action of a. The term "BPS 
equation" derives from the observation that (2.45) inserted in (2.23) implies the saturation of the 
BPS bound (2.24), which simplifies to 

E[P] = 8ttQ[P] = 87rTr(Pa(l-P)a t P) = 8tt Tr (P a a f P - a P a)) . (2.48) 

For BPS solutions of the particular form (2.33) we indeed recover that E[P] = 8irTrP'. Clearly, 
these BPS solutions constitute the absolute minima of the energy functional within each Grass- 
mannian. 

When the parametrization (2.36) is used, the BPS condition (2.45) simplifies to 

a\T-) = \Tf}j/ for some rxr matrix 7 = (7/) (2.49) 

which represents the action of a in the basis chosen for imP. For instance, the ket \Tq) in (2.43) 
indeed obeys 



a\T Q ) = \T Q ) 1Q with 7 Q 



( VI ... \ 

V2 



(2.50) 



• ... VQ^T i 

\o ... / 



and the diagonal projector P^ in (2.42) is BPS, but it is by far not the only one. 

Any basis change in imP induces a similarity transformation 7 1— > ^r -1 , which leaves P 
unaltered and thus has no effect on the value of the energy. Therefore, it suffices to consider 7 to 
be of Jordan normal form. More generally, a unitary transformation (2.41) is compatible with the 
BPS condition (2.49) only if 

U^aU\T) = \T)~fu for some rxr matrix 7c/ , (2-51) 

which implies that 

uj\T) = \T) luj with lu) = lu - 1 . (2.52) 

The trivially compatible transformations are those in U(imP) xU(kerP), which can be subsumed 
in S G GL(imP) and lead to r yjj=T~ lr )T . Another obvious choice are rigid symmetries of the 
energy functional, as given in (2.13), (2.16), and (2.11) for W=V* . The challenging task then is to 
identify the nontrivial BPS-compatible unitary transformations, since these relate different a-stable 
subspaces of fixed dimension in H® n and thus generate the multi-soliton moduli space. For \Tq) 
in (2.43) we shall accomplish this infinitesimally in Section 3.3. 

2.4.1 Abelian solitons 

Let us take a closer look at the BPS solutions of the noncommutative U(l) sigma model. All 
finite-energy configurations are based on P=0 and have rank r = Q < 00, thus E = 8nr. The task 
is to solve the "eigenvalue equation" 

\ 



with a s G C , (2.53) 



a\T) = |T> 7 for 7 = 



'a s 1 
a s 1 



s =l \ a 3 
... 
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where the Jordan cells have sizes r s for s = 1, . . . , q, with Y11=i r s = r - For a given rank r, the above 
matrices 7 parametrize the r-soliton moduli space. The general solution (unique up to cell-wise 
normalization and basis changes \T^) — > |T( s ))r( s )) reads 

|T) = (|TW> ... |TW>) with |T«> = (|a s ) a t|a.) ... ^(atj^V,)) (2-54) 

and is based on the coherent states (2.14). In the star-product picture, the corresponding $ repre- 
sents r lumps centered at positions a s with degeneracies r s in the xy plane. Lifting a degeneracy 
by "point-splitting", the related Jordan cell dissolves into different eigenvalues. Hence, the generic 
situation has r s =l Vs, and so 

7 = diag(ai,a 2 ,...,a r ) |T) = (jai) |a 2 ) ... K)) , (2.55) 

which yields the projector 

^ = E !«*>(<«») <«*l • (2-56) 

Each solution can be translated via a unitary transformation mediated by D([3), which shifts 
an 1 ^ a^+/3 W, and rotated by which moves a.£ i-> e l1? c^ W. The individual values of (the 

soliton locations) may also be moved around by appropriately chosen unitary transformations, so 
that any r-soliton configuration can be reached from the diagonal one, which describes r solitons 
on top of each other at the coordinate origin [9]: 



r-1 



\T) = U\T r )Y P = UP r U^ with P r = ^ \m)(m\ . (2.57) 

m=0 

We illustrate the latter point with the example of r = 2. Generically, 

(|ai) |<*2» = US (\0)\1)) = U(\0)\1))T with US = |c*i)(0| + |a 2 )(l| + . . . , (2.58) 
where the omitted terms annihilate |0) and |1). Factorizing US yields 

1 /. . . ,\ fe~^ \ /-sin/3' cos /A /(0|\ 

and T = 

with a = (ai|a 2 ) = e " ia2 -5l«il 2 -5l^l 2 = e^ 7 " 7 ') cos(/?-/3') , 
so that \a\ 2 = e - !" 1- " 2 ! 2 , and /3+j3' and 7+7' remain undetermined. The rank-2 projector becomes 

P = z^r^(\a 1 ){a 1 \ + \a2){a2\-a\a l ){a 2 \-a\a 2 ){a 1 \) = U (|0)<0| + |1)(1|) U* . (2.60) 



cos/3 cos /A /e i7 ^ (2.59) 
sin/? sin/3' J I e i7 ' 



2.4.2 Nonabelian solitons 



We proceed to infinite-rank projectors. For simplicity, let us discuss the case of U(2) solitons - 
the results will easily generalize to U(n). Clearly, we can imbed two finite-rank BPS solutions 
(with R=0) into \]{H®H) by letting each act on a different copy of H. Such configurations are 
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noncommutative deformations of the trivial projector P = ((jo) ano - thus represent a combination of 
abelian solitons. Therefore, we turn to projectors with R = 1 and so, formally, r = oo + Q. Given 
(2.33) such solutions may be considered as noncommutative deformations of the U(2) projector 
P = (oo)- F° r t ms reason, one expects the generic solution |T) to (2.49) to be combined from a 
full set of states for the first copy of TC and a finite set of coherent states \ag) in the second copy 
of H, 



|T) 



\H) 







\ai) \a 2 ) ... |«q) 



kl=l 



with P = 87r<5. This projector is of the special form (2.33), with P = (JJJ) and U = 1. By a 
unitary transformation on the second copy of TL such a configuration can be mapped to the diagonal 
form 

(\H) ... \ 



(2.62) 





JO) |1) 



It is convenient to reorder the basis of imP^ such that 

|r d > = 

where 



\H) 

IQ-i) , 



m 



(2.63) 



(2.64) 



m=0 



m=Q 



denotes the Qth power of the shift operator. The form of (2.63) suggests to pass from states 
\T) = (j^j) with R =1 to operators f = 0^ on H: 



\T) = f\H) => P = f (ftfj-ift 
In fact, it is always possible to introduce T as 



(2.65) 



£ = £|t/)(£-i| 



Ti\£-1) for t = l,2 and £=l,...,r. (2.66) 



We may even put T^T = In by using the freedom T — ► f f with an operator f = (T^T) 1 / 2 . Our 
example of |Td) in (2.63) is already normalized since 



S 1 S\ = l n but 5} Si = 1«-|0)<0| 



SqSq + Pq = In • (2.67) 



It is instructive to turn on a BPS-compatible unitary transformation in (2.33). In our exam- 
ple (2.63), we apply [25] 



Ufa) 



(s q\J s h s q^nq+„h 



s. 



t fj 



'Nr. 



\ y/NQ+flfi Q y/NQ+flji J 



with N Q = a) Q a Q = N(N-l) • • • (N-Q+l) 

(2.68) 
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to f d = (5q,Pq)* of (2.63). With the help of 

S Q P Q = = N Q P Q and Sq^N^ = a Q (2.69) 

we arrive at 

*M = W ft = (; S ) = (;°) f _: f(,)f . (2.70, 

This transformation can be regarded as a regularization of since 

limf(^) = il%)f A with <5=limarg/x and thus MmPUi) = P d , (2.71) 

fx— »0 u /i— >0 ^-^0 

but note the singular normalization in the limit! For completeness we also display the transformed 
projector, 



(2.72) 



How do we see that such projectors are BPS? Writing \T) = T\Tt), the BPS condition 

a\lf) = \Tf)j/ implies a% = ftf with i = l,2 (2.73) 

for some operator 7 in We do not know the general solution for arbitrary 7. However, an 
important class of solutions arises for the choice 7 = where the BPS condition reduces to the 
"holomorphicity condition" 9 

[a,fi]=0 for i = 1,2 . (2.74) 

These equations are satisfied by any set of functions {Ti,^} of a alone, i.e. not depending on a'. 
Indeed, for the example of in (2.63), we concretely have 

7 = aP Q + (l w -p g ) a ^2, (2.75) 

while T(/x) in (2.70) is obviously holomorphic and thus BPS. Because T(/x) emerges from via a 
BPS-compatible unitary transformation it shares the topological charge Q and the energy E = 8irQ 
with the latter. 

The generalization to arbitrary values of n and R < n is straightforward: T becomes annxi? 
array of operators (T^), on which left multiplication by a amounts to right multiplication by an 
RxR array of operators 7. For the special choice 7^ = 5^a, any collection of holomorphic 
functions of a serves as a solution for T^. Quite generally, one can show that for polynomial 
functions T t L (a) the resulting projector has a finite topological charge Q given by the degree of 
the highest polynomial and is thus of finite energy [26, 5, 27]. In this formulation it becomes 
evident that nonabelian solutions have a smooth commutative limit, where y/29a — > z and — > 0. 
Indeed, they are seen as deformations of the well known solitons in the Gr(n, R) = u(R)xU(n-R) 
Grassmannian sigma model. 10 Hence, the moduli space of the nonabelian solitons coincides with 
that of their commutative cousins. By rescaling T — > TT with a C-valued RxR matrix T we 

9 This may even be the general case: If there exists an invertible operator T solving aF — Tj, then the general 
solution to (2.73) reads Ti = TJ with [a, Ti] = 0, and T can be scaled to unity. 
10 For the above-discussed example one has Gr(2, 1) = CP 1 . 
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can eliminate R complex parameters from nR independent polynomials. For a charge-Q solution, 
there remain nRQ + (n— 1)R complex moduli, of which (n— 1)R parametrize the vacuum and nRQ 
describe the position and shape of the multi-soliton [26]. For the case of CP 1 this yields a complex 
2Q dimensional soliton moduli space represented by = ( a "q+"+^)- 0m* sample calculation 

above suggests that taking T2 — ► fj, and then performing the limit \i — > one recovers the complex 
Q dimensional moduli space of the abelian solitons given by \T) = (|«i) . . . |ttQ_i)) as a boundary. 

2.5 Some non-BPS solutions 

For the record we also present a particular class of non-BPS (and non-anti-BPS) Grassmannian 
solutions to the equation of motion [AP, P] = 0. From the action of A on a basis operator, 

A\m)(n\ = (m+n+1) \m)(n\ — \Jmn \m— l)(n-l| - V(m+l)(n+l) \m+l){n+l\ , (2.76) 

we infer that A maps the A;-th off-diagonal into itself, so in particular it retains the diagonal: 

A\m)(m\ = (2m+l) |m)(m| — m \m— l)(m— 1| — (m+1) |m+l)(m-|-l| . (2.77) 

It follows that every diagonal projector is a solution. Let us first consider the case of U(l). Given 
the natural ordering of the basis {\m)} of H, any diagonal projector P' of finite rank r=Q can be 
written as 

q r s -l q 

P' = ^2^2\m s +k)(m s +k\ = ^ P Te S ms with m s+1 > m s +r s Vs (2.78) 
s=l k=0 s=l 

and Ylt=i r s = r - Via 

AP' = m s \m s )(m s \ - m s \m s -l)(m s -l\ 

7^1 1 (2.79) 

- (m s +r s ) \m s +r s ){m s +r s \ + (m s +r s ) |m s +r s -l)(m s +r s -l|| 
its energy is easily calculated as 

q 

±E[P'\ = Q + m s , (2.80) 

s=l 

which is obviously minimized for the BPS case q = 1 and mi = 0. The energy is additive as long 
as the two projectors to be combined are not getting "too close". This picture generalizes to the 
nonabelian case by formally allowing r s and m s to become infinite. In particular, the energy does 
not change when one imbeds P' (or several copies of it) into \J(Ti.® n ) and adds to it a constant 
projector as in (2.33). Hence, with the proper redefinition of the m s , (2.80) holds for nonabelian 
diagonal solutions as well. 

The inversion P — > 1— P generates additional solutions, which for the structure in (2.33) and 
BPS-compatible unitaries U are represented as 

P = U(P-P')U^ with P= (l n -P)®l w and (l-P)P' = = P'(l-P) , (2.81) 

When t—P+P' is BPS, then P becomes anti-BPS with topological charge Q[P] = — TrP', produc- 
ing an anti-soliton with energy E[P\ = 87rTrP'. It is possible to combine solitons and anti-solitons 
to a non-BPS solution via 

P = Psol + P^i Provided P sol P- { = = P^P sol , (2.82) 
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so that their topological charges and energies simply add to 

Qi p ] = Q[P S oi] + Ql P ^] aild E[P] = E[P sol ]+E[P^]. (2.83) 

For the diagonal case, this is included in the solutions discussed above. Examples for U(l) and 
U(2) (which appeared in [28]) are (for m>r) 

m—l 

P = In-Pra + Pr = 1<H ~ £ \ k )( k \ 

k=r 

P = P ri ® (Iw-Pra) 

respectively. Besides the global translations and rotations, other unitary transformations are con- 
ceivably compatible with the equation of motion [AP, P] = 0, generating moduli spaces of non-BPS 
Grassmannian solutions. Outside the Grassmannian manifolds, many more classical configurations 
are to be found. 

3 Fluctuation analysis 

In order to investigate the stability of the classical configurations constructed in the previous section, 
we must study the energy functional E in the neighborhood of the solution under consideration. 
Since the latter is a minimum or a saddle point of E, all the linear stability information is provided 
by the Hessian, i.e. the second variation of E evaluated at the solution. The Hessian is viewed 
as a linear map on the solution's tangent space of fluctuations, and its spectrum encodes the 
invariant information: Zero modes belong to field directions of marginal stability (and extend to 
moduli if they remain zero to higher orders) while negative eigenvalues signal instabilities. A 
perturbation in such a field direction provides (part of) the initial conditions of a runaway solution 
in a time extension of the model. We cannot be very specific about the stability of a general 
classical configuration. Therefore, we shall restrict our attention to the stability of BPS solutions 
(as introduced above), at which the Hessian simplifies sufficiently to obtain concrete results. Since 
any BPS configuration is part of a moduli space which is imbedded in some Grassmannian which 
itself lies inside the full configuration space of the noncommutative U(n) sigma model, the total 
fluctuation space contains the subspace of Grassmannian fluctuations which in turn includes the 
subspacc of BPS perturbations, the latter being zero modes associated with moduli. In order 
to simplify the problem of diagonalizing the Hessian we shall search for decompositions of the 
fluctuation space into subspaces which are invariant under the action of the Hessian. As we consider 
configurations of finite energy only, admissible fluctuations (f) must render S 2 E finite and keep the 
"background" <I> unitary. Furthermore, they need to be subject to the same conditions as $ itself: 
<fi is bounded, [a,(f>] and [at,<f>] are Hilbert-Schmidt, and A</> is traceclass. Finally, in keeping with 
our restricted notion of Grassmannian, we do not admit hermitian fluctuations 11 which alter the 
topological charge. 

3.1 The Hessian 

The Taylor expansion of the energy functional around some finite-energy configuration $ reads 

= E ^ + /*« £$1*1 *M + I M^'^wtBn?) 1 * 1 *W +•■■ (31) 

=: g[£] + + £ (2 ) + ••• , 

11 Perturbations inside the Grassmannian must be hermitian (see Section 3.2 below). 



with Q = r—m and J-E = r+m , 

^ 8 " (2.84) 

with Q = r\—T2 and ^-E = r\+r2 , 
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where the U(n) traces are included in J d 2 z. The perturbation <j> is to be constrained as to keep 
the background $ unitary. Since we compute to second order in <f>, it does not suffice to take 
4> £ X$U(W ffin ). Rather, we must include the leading correction stemming from the exponential 
map onto U(7Y® n ), which generates the finite perturbation cf> = <£>. The latter is subject to the 
constraint 

($ t +0 t )($+<£) = 1 =► </»t = _$t^$t + $t^$t^$t + O(0 3 ), (3.2) 

which we shall use to eliminate 0^ from the variations. It is important to realize that in this way the 
term linear in $ generates a contribution to E^[$>, <j)\. Performing the expansion for the concrete 
expression (2.8) and using (3.2) we arrive at 12 

EW [&,</>] = TrTrjfat, $V$ f ] + [a\ $] [a, $V$ f ] - [a f , $*] [a , 0] - [a f , 0] [a, 

= 2vrTr{(A$ t $- $ f A$) $ f 0} = 0, (3.3) 

£( 2 )[$,^] = TrTrlla 1 ", $V$ f ] [a, 0] + [a f , <j>\ [a, $V$ f ] 

- [a f , $V$t0$t] [a, $] - [a f , $] [a, $V$V$t]} 

= 2vrTr{$ t (/>$ t </>$ t A$ - &<j>& A<^>} 

= 2^Tr{0 t A0-0 t ($A$ t ) 0} + O(^ 3 ) =: 2tt Tr{<ft H <j)} + O(0 3 ) , (3.4) 

defining the Hessian H = A — ($A3>t) as a self-adjoint operator. Hence, our task is essentially 
reduced to working out the spectrum of the Hessian. Since A is clearly a positive semidefinite 
operator, an instability can only occur in directions for which (<&A<&t) is sufficiently large. 

For later reference, we present the action of H in the oscillator basis, 
H ^2^,e\m}(£\ = J2( H ^m,e\m)(£\ with 

m,£ m,£ 

(H(f>) mt e = {m+i+l) (j) m/ - v 7 (m+l)(£+l) (j) m +i,e+i - Vm£^ m -i,e-i ( 3 - 5 ) 

" {(j + k+l)$ jtk -y/(j+l)(k+l)$ j+1>k+1 - v9fc*j-i,fc-i}^M > 

j,k 

where Q m ,£ as well as <\> m ,i are still nxn matrix- valued. At diagonal abelian backgrounds, as given 
in (2.78), the matrices reduce to numbers and the latter expression simplifies to 

{H(f>) mie = (m+£+l - 2b m ) (j) m/ - y/(m+l)(e+l) <f>m+i,e+i — Vrnl ' <f> m -i t g-i with 

i (3.6) 

(m s +r s ) + 

$m,m s +r s — 1 

(m s +r s )} , 

s=l 

which for fixed t differs from (A<^>) m ^ in at most Aq entries. For diagonal U(l) BPS backgrounds 
$ = cf> r = 1 — 2P r this reduces further to 

b m = r (S m: r-1 + S m ,r) ■ (3.7) 

It is important to note that these expressions do not yet yield a matrix representation of H because 
the constraint (3.2) on the allowable perturbations must still be taken into account. We can do 
this by replacing <p with <j> — Q^ffr everywhere; then (f> becomes unconstrained. 

12 The same result is obtained by considering E[exp{t(f>}&] up to 0(t 2 ). 
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3.2 Decomposition into even and odd fluctuations 

We specialize to Grassmannian backgrounds, <E> = 1 — 2P = <J>t, characterized by a hermitian 
projector P and obeying $ 2 = 1. Any such projector induces an orthogonal decomposition 

C n ®H = P(<D n ®H)®{l-P)(*D n ®H) =: imP © kerP , (3.8) 

and a fluctuation </> decomposes accordingly as 

(j) = P(j)P + (l-P)0(l-P) + P(f>(l-P) + (l-P)^P, (3.9) 

V v ' V v ' 

where the subscripts refer to "even" and "odd", respectively. Since <I> acts as —1 on imP but as 
+ 1 on kerP, we infer that 

®<fi c = </> e <£> and $ O = -00$ => <t>l = -<f> e and 4>\ = <f> (3.10) 

to leading order from (3.2), i.e. even fluctuations are anti-hermitian while odd ones are hermitian. 
This implies that odd fluctuations keep $ inside its Grassmannian, but even ones perturb away 
from it. It also follows that in 

Tr(0t A<f>) = Tr(4 Acb c ) + Tr(4 A0 O ) + Tr(4 A0 O ) + Tr A0 C ) (3.11) 

the last two terms cancel each other. Furthermore, the equation of motion (2.44), [AP, P] = 0, 
implies that 

($A$) = Tr (<j>l $A$ <f> ) = = Tr(0 o $A$ (f> e ) , (3.12) 

because only an even number of odd terms in a product survives under the trace. Combining (3.11) 
and (3.12) we conclude that 

[<J>,0c+0o] = P (2) [$,0c]+P (2) [$,0o] , (3.13) 

which allows us to treat these two types of fluctuations separately. 

The above decomposition has another perspective. Recall that any background configuration 
$ being unitary can be diagonalized by some unitary transformation, 

$ = ?7$ d [/t = [/diag({e iAl }) U ] . (3.14) 

When $ is hermitian, i.e. inside some Grassmannian, the diagonal phase factors can be just +1 
or —1, and U is determined only up to a factor V € U(imP) x U(kerP) which keeps the two 
eigenspaces imP and kerP invariant. Adding a perturbation <j> lifts the high degeneracy of so 
that the diagonalization of <I>+0 requires an infinitesimal "rotation" K of imP and kerP inside 7i 
as well as a "large" rediagonalization V inside the two eigenspaces. Modulo higher order terms we 
may write 

$ + = U(l+K)V($ d + <t> d )V*(l-K)U* = $ + U{V(f> d V^ + [K , $ d ]}tf 

(3.15) 

with [V, <J>d] = and K = — infinitesimal , 

where 4> d is a purely diagonal and anti-hermitian fluctuation. Since V depends on <f> it should not 
be absorbed into U. It rather generates all non-diagonal fluctuations inside imP and kerP, allowing 
us to rewrite 

VfaV* = <f/ d + [A e ,4>d] , with At = -A e (3.16) 
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being a generator of U(imP) x U(kerP) and a modified diagonal perturbation c/)' d . Redenoting also 
K = eA Q with a real and infinitesimal e and a generator A G of the Grassmannian, the general 
fluctuation is parametrized as 

<f> = U{<f>' A + [A c ,^ d ] + [A , e$ d ]}C/t (3.17) 

and decomposed (after diagonalizing the background via U) into a "radial" part 4>' d and an "angular" 
part <p & = [A, any] with A generating \J(7i® n ) [29]. For a Grassmannian background all terms have 
definite hermiticity properties, and we can identify 

<Pe = U{<p' d + [Ae,^]}^ and <p o = U{[A Q , e$ d ]}rf = e [UA tf , $] . (3.18) 

We have seen in (3.13) above that the even and odd fluctuations can be disentangled in E^ 2 \ It 
is not clear, however, whether the diagonal perturbations can in turn be separated from the even 
angular ones in the fluctuation analysis (but see below for diagonal backgrounds where U = 1). 

3.3 Odd or Grassmannian perturbations 

As far as stability of BPS configurations is concerned, the odd perturbations are easily dealt with 
by a general argument. Since a shift by <f> a keeps <3> inside its Grassmannian, wherein already 
minimizes the energy, such a perturbation cannot lower the energy any further and we can be sure 
that negative modes are absent here. Therefore, solitons in the noncommutative Grassmannian 
sigma model are stable, up to possible zero modes. An obvious zero mode is generated by the 
translational and rotational symmetry. A glance at (2.13) and (2.16) shows that the corresponding 
infinitesimal operators A are given by 

Atrans = a — a a and A rot = id a) a , (3.19) 

which indeed leads to the annihilation of [A trans , <E>] and [A rot , <£] by H as we shall see. 

From the discussion at the end of Section 2.4 we know that for r>l there are additional zero 
modes inside the Grassmannian, because the multi-soliton moduli spaces are higher-dimensional. 
Parametrizing these BPS-compatible perturbation as follows, 

$ + o = E/ B $E4 = $ + e[A B ,$] + 0(e 2 ) for U B = e eA B with A B = -A B , (3.20) 
the corresponding Lie-algebra element (2.26) becomes 

u = ul[a,U B ] = ( e — dA B-l)a = e [a , A B ] + 0(e 2 ) , (3.21) 
and the condition (2.52) of BPS compatibility to leading order in e reads 

[a,A B ]|T) = |T)7 aJ for some rxr matrix 7^ . (3.22) 

Let us try to find A B in the abelian case by perturbing around the diagonal BPS configuration 

r-l 

\T r ) = (|0> |1) ... |r-l» $ r = l n ~ 2^|m)H . (3.23) 

m=0 

Expanding the generator 

A B = J2A m/ \m){£\ with A m/ = -A £im (3.24) 

m,£ 
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it is easy to see that the even part of A B automatically fulfils (3.22), and so restrictions to A m ^ 
arise only for the odd components. In fact, only the terms with m>r and £<r in the above sum 
may violate (3.22), which leads to the conditions 

%/m+l A m+ i t £ — VI A m< g_i = for all m> r and I < r . (3.25) 

Since £=0 yields A m+ i ) o=0 as a boundary condition, this hierarchy of equations puts most compo- 
nents to zero, except for m = r, . . . ,r+£ at any fixed £<r. These remaining r(r+l)/2 components 
are subject to r(r— 1)/2 equations from (3.25), whose solution 

W+j = ^4)trV 3 ~Xt^) K ^ fOT j = l>-,r-l-t at t< r ~ 2 ( 3 - 26 ) 

fixes r{r— 1)/2 components in terms of the r— 1 components appearing on the right hand side, 
which therefore are free complex parameters. The rth free parameter A r r _i does not enter and 
is associated with the rigid translation mode. 13 Finally, the ensueing BPS perturbation (3.20) is 
found to be 

oo r— 1 

<f> = e[A B ,<D r ] = e^^(A m ,,|m)(^| + A m ^)(m|) , (3.27) 

m=r £=0 

with A m £ taken from (3.26). To higher orders in e, the BPS-compatibility condition is not automat- 
ically satisfied by our solution (3.26) but this can be repaired by adding suitable even components 
to A B . Our result ties in nicely with the observation of r complex moduli in (2.56) whose shifts 
produce precisely r complex zero modes. For the simplest non-trivial case of r=2, one can also 
extract these modes from differentiating (2.59) with respect to a.\ or oli- 

3.4 Even or non-Grassmannian perturbations 

The stability analysis of BPS configurations inside the full noncommutative U(n) sigma model 
requires the investigation of the even fluctuations 4> c as well. Here, we have only partial results to 
offer. 14 Yet, there is the following general argument which produces an unstable even fluctuation 
mode neg (but not an eigenmode) for any noncommutative multi-soliton $ = 1 — 2P with Q > 0. 
For this, consider some other multi-soliton <3> = 1 — 2P which is contained in $ in the sense that 

im(P) C im(P) PP = PP = P . (3.28) 

We then simply say that P C P. It follows that their difference II is the orthogonal complement 
of P in im(P), 

U = P - P C P => n 2 = II and nP = = PII. (3.29) 
In particular, we may choose P = 0. For any such pair (P, P) there exists a continuous path 
$(a) = e isn (l-2P) = 1 - 2P + (l-e is )n = 1 - 2P - (l+e is )n 

(3.30) 

connecting $(0) = $ = 1-2P with $(7r) = $ = 1-2P . 

Note that 3>(s) interpolates between two different Grassmannians, touching them only at s=0 and 
s=ir. Since we assumed that and <I> are BPS we know that 

E[§] = 8irQ and E[$] = 8ir Q (3.31) 

13 The rigid rotation mode is absent because $ r is spherically symmetric. 
14 Even for the commutative sigma model this is an open problem [2]. 
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with Q and Q being the topological charges of P and P, respectively. 15 Inserting (3.30) into the 
expression (2.8) for the energy and abbreviating 1— e ls = p we compute 

E(s) := E[$(s)] = 8^Q-2vr( / 9+p)Tr([a,P][n,a t ] + [a t ,P][n,a]) +2vr / 9pTr([a,n][n,a t ]) 

= 8vrQ-47r(l-coss)Tr([a,P][n,a t ] + [a^P][U,a] - [a,U}[U,a^]) , (3.32) 

where the two traces could be combined due to the relation 

p + p = pp = 2(1— cos s) . (3.33) 

Luckily, we do not need to evaluate the traces above. Knowing that E(tt) = 8irQ we infer that the 
last trace in (3.32) must be equal to Q—Q and hence 

E(s) = 8vr Q - 4tt(1- cos s)(Q - Q) = 8tt(Q cos 2 f + Q sin 2 f ) 

(3.34) 

= 8ttQ - 2tt{Q -Q)s 2 + 0(s 4 ) . 

Evidently, for Q > we can always lower the energy of a given BPS configuration by applying an 
even perturbation cp ncg = — iell towards a BPS solution with smaller charge Q < Q. The higher 
the charge of P the more such modes are present. Yet, they are not independent of one another 
but rather span a cone extending from the background, as we shall see in examples below. In fact, 
there is no reason to expect any of these unstable fluctuations to represent an eigenmode of the 
Hessian, and in general they do not. Nevertheless, their occurrence again demonstrates that there 
must be (at least) one negative eigenvalue of H, and for diagonal U(l) BPS backgrounds we shall 
prove in Section 4.3 that there is exactly one. The only stable BPS solutions are therefore the 
"vacua" defined by P'=0 in (2.33) and based on a constant U(n) projector P. This leaves no stable 
solutions in the abelian case besides $ = 1^. 

In addition to the unstable mode, there exist also a number of non-Grassmannian zero modes 
around each BPS configuration, which generate nearby non-BPS solutions to the equation of mo- 
tion. This will become explicit in the examples discussed below. 

4 Perturbations of U(l) backgrounds 

4.1 Invariant subspaces 

We specialize further to diagonal U(l) backgrounds $ as given by (2.78) (not necessarily BPS). 
Using (3.5) it is easy to see that H maps any off-diagonal into itself. Let us parametrize the k-th. 
upper diagonal as 

oo 

h+k) = M(fe)m \m)(m+k\ with p {k)m e <D for k = 0, 1,2, ... . (4.1) 

m=0 

The hermiticity properties (3.10) of the perturbations demand that we combine T>k an d into a 
subspace £ & of the fc-th upper plus lower diagonals by defining 

£k ■= {<t>(k) | (t>(k) = 4>(+k) ~ ^4>\ +k) ^} , which implies 4>\ k) = -$0 (fc) $. (4.2) 

The direct sum of all is the full admissible tangent space to U(H), and so we may decompose 

oo oo oo 

= = ^2^2{^(k)m \m)(m+k\ =F U(k)m \m+k)(m\} , (4.3) 

k=0 t=0m =0 

15 Note that Q need not be smaller than Q when r(P) is infinite. 
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with the sign depending on whether the component is even or odd. Since H maps £ k into itself, 
the bilinear form defined in (3.4) is block-diagonal on the set of 

2vrTr{0{ fe) ^0 (£) } ~ 5 M , (4.4) 

which implies the factorization 

E {2) [$,E k <t>(k)} = E fc ^ (2) [*^(fc)] • (4-5) 

In other words, forms an .ff-invariant subspace for each value of k. In particular, £$ is the space 
of admissible skew-hermitian diagonal matrices, and the (purely imaginary) diagonal fluctuations 
0(0) = (j>d G £o can be considered on their own. 

4.2 Results for diagonal U(l) backgrounds 

Apparently, for diagonal U(l) backgrounds it suffices to study the second variation form E^l®, .] 
on each subspace £ k (k > 0) separately. Let us give more explicit formulae for the restriction 
of E^ [<£,.] to these subspaces. To this end, we denote the non-zero entries of <3? by 8j := $jj 
(j = 0, 1, 2, . . . ). We have 5j = ±1 for each j, and the set 

J := {j > | S j+1 + Sj} (4.6) 

is finite. Then, a straightforward calculation shows that the restriction of E^[<&, .] to 



£o = {4>(o) | 0(0) = i^2<t>m\m)(m\ with (j) m £ R} is given by 

m=0 

oo 

±eW[<S>^ (0) ] = Y J (rn+l)(<t>m + i-0m) 2 -2j2(j+l)(<P] + i+4> 2 3 ) ■ 



(4.7) 



The formula for E^ [<£, .] on £k>o is more complicated. Considering a fixed k-th upper diagonal 
Pfc as parametrized in (4.1), we first evaluate (suppressing the subscript (k)) 

oo oo 

Tr{</>{ +fc) A(f) {+k) } = Wm+lfi m+1 - Vm+k+ln m \ 2 = \k \^\ 2 + \ ^ Rm(n) , (4.8) 

m=0 m=0 

with 

^m(M) : = |Vm+l/x m+ i - Vm+k+1 Li m \ 2 + | Vm+fc+l fi m +i ~ Vm+ln m \ 2 (4.9) 
= (Vm+k+l - v / m+l) 2 (|/x m+ i| 2 + |/U m | 2 ) + 2\/ m+k+l\J m+1 |/x m+ i - ;U m | 2 . 

Armed with this expressions, we compute the second variation on To state the outcome, it is 
useful to introduce the index sets 

J-k := {j-ke N |j € J} and A := (J U (J - fc)) \ (J n (J - fc)) , (4.10) 

i.e. ^4 is the symmetric difference of J and J—k. For </>( fc ) G a direct calculation results in 

±EW[$,<p {k) ] = k\^\ 2 + E ^-oo + E^'+^d^i' + i^i 2 ) ( 4 - n ) 

jeN \A jeA 



2^(j+l)(|^| 2 + |^ +1 | 2 ) - 2 £ (j+k+l)\ H \ 2 - 2 £ 0"+fc)lA*. 
jeJ jeJ-fc jeJ-fc+i 



20 



This expression simplifies when k is greater than the largest element of J+l. Then the sets 
J—k and J— k+1 are empty, so A = J, and (4.11) takes the form 



+ 



+ ^E(n 2 + iwi 2 ) • ( 4 - 12 ) 

One sees that for each je J the corresponding coefficient Hj decouples from all successive coeffi- 
cients /U m >j. Since the elements of J signify in the string of fi m the boundaries between even and 
odd fluctuations, this observation confirms the decomposition (3.13) of E^ into an even and odd 
part also after the restriction to £fc. 16 Furthermore, we note that Rj(n) > unless fij = Hj+i = 0, 
which makes it obvious from the expression (4.12) that the quadratic form E^ is strictly positive 
on each E k with k > maxjgj^'+l). For the remaining £k we have to work a little harder. 

4.3 Results for diagonal U(l) BPS backgrounds 

The idea is to pursue the reduction of E^ to a sum of squares whenever possible. For diagonal 
BPS solutions (2.57) 

<$> r = l n - 2V|m)(m| = V<S m |m)H with 5 m = { , (4.13) 

^— 1 ^— ' +1 for m > r 

m=0 m=0 K — 

this strategy turns out to be successful at all k > 1 but breaks down at k = 0. We note that now 
J = {r— 1}, which implies a distinction of cases: "very off-diagonal" perturbations have k > r, 
"slightly off-diagonal" perturbations occur for 1 < k < r, and diagonal perturbations mean k = 0, 
to be discussed last. It will also be instructive to visualize the Hessian on the space £ k in matrix 
form, 



£ (2) [$r,0( fc) ] ^ 

defining an infinite-dimensional matrix 



27T ^ VmH^M 

m,e=o 



(4.14) 



4.3.1 Very off-diagonal perturbations 

Since J consists just of one element, at k > r the matrix H ^ splits into two parts only, of which 
the odd one has finite size r. More explicitly, 



H (k) 



H, 



(fc) 

Gr(P) 



with 



/ k + 1 



1 rr(k) 
2- f2 Gr(P) 



-VW+ij 

y/l{k+l) k + 3 -y/2(k+2) 
-y / 2(k+2) fc + 5 



(k) 

kerP ' 



(4.15) 



2r + fe-3 -y/( r -l)(r+k-l) 



-y/(r-l)(r+k-l) 2r + k-l-r j 

(4.16) 



'For | J| > 1 the even and/or odd part of is split further. In total, decomposes into at least | J\ + l blocks. 
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1 n 

2 rL kcrP 



(k) _ 



( 2r + k + l-r -^(r+l)(r+k+l) 

-y/(r+l)(r+k+l) 2r + fc + 3 -y/(r+2)(r+k+2) 

-^/(r+2)(r+k+2) 2r + k + 5 

V 



(4.17) 



where the boldface contributions disturb the systematics and originate from the <!>A<E>t term in the 
Hessian. The previous paragraph asserts strict positivity of H for k > r. Due to the finiteness of 
^Gr(P)' Hessian on 6k> r features precisely r positive eigenvalues 17 in its spectrum. We shall 

(k) 

argue below that H^ p contributes a purely continuous spectrum R + for any k. Thus, the above 
eigenvalues are not isolated but imbedded in the continuum. 

It is instructive to look at the edge of the continuum. The non-normalizable zero mode of 

(k) 

^Gr(P) * s explicitly given by 



zero 



(r+k+l)(r+k+2) ■ ■ ■ (r+k+m) 
(r+V)(r+2) • • • (r+m) 



(r+m+l)(r+m+2) ■ ■ ■ (r+m+k) 
Vr\l — rTT77~T7T\ ~ , r,\ > \ AA °) 



(r+l)(r+2)-- - (r+k) 



which grows like m k / 2 when m = 0,1,2,... gets large. Being unbounded for &:>0 this infinite 
vector does not yield an admissible perturbation of <3> r however: 5 2 E is infinite. 

4.3.2 Slightly off-diagonal perturbations 

For each value of k in the range 1 < k < r we already established in Section 3.3 the existence of an 
odd complex normalizable zero mode connected with a moduli parameter. Nevertheless, as we will 
show now, E^[& r , 4>{k)\ remains positive for 1 < k < r albeit not strictly so. In order to simplify 
(4.11) we make use of the following property for Rj(fi) as defined in (4.9): 



l-i 



l-i 



^Rjifi) = k\m\ 2 - k\/j, m \ 2 + 2 ^ Wj+ 1 H+i ~ Vj+ k + 1 H\ 



j=m 



l-l 



(4.19) 



k\fi m \ 2 - k\fii\ 2 + 2 ^ W Hj+i ~ 



]=m 



Here, by definition, all integers are non-negative and all sums are taken to vanish if m > /. The 
two equations above are easily proved by induction over I, starting from the trivial case I = m. We 
now employ the algebraic identities (4.19) to rewrite (4.11) for 1 < k < r and obtain 



±EW[$ r ,<f> {k) ] = k\^\ 2 + X>» 



(4.20) 



r-k-2 r-2 
j=0 j=r-k 



17 meaning that the corresponding modes are normalizable eigenvectors, i.e. Hilbert-Schmidt 
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which is indeed positive semi-definite. We observe that the even part of now consists of two 
disjoint pieces, containing {fig, . . . ,/x r -fc-i} and {fij> r }, which are separated by the odd perturba- 
tions {^ r _fc, . . . ,/J, r -i}. If the right-hand side of (4.20) is equal to zero, then all components /j,j> r 
must vanish because each Rj is a strictly positive quadratic form of jij and fJ,j+i- In contrast, the 
remaining coefficients /Uj <r need not vanish at the zero set of 4>ha[. In order to find the 

zero modes of the Hessian, it is convenient to visualize it again in matrix form, namely 

where the blocks have sizes r—k, k, and oo, respectively. From (4.20) we extract 



2 I1 imP 



(fc) _ 



/ fc+1 -y/l(k+l) 

-y/l(k+Tj k + 3 - v / 2(A;+2) 
- v / 2(fc+2) k + 5 



2r-k-3 -y/( r -k-l)(r-l) 

-y/(r-k-l)(r-l) 2r-k-l-r J 

(4.22) 



1 rr(fc) 
2- f2 Gr(P) 



/ 2r - fc + 1 



V(r-fc+l)(r+l) 



-V(r-A;+l)(r+l) 2r-fe + 3 



2r + /c-3 —y(r-l)(r+A;-l) 



V(r-l)(r+fc-l) 2r + A;-l-r / 

(4.23) 



and H^2 P being identical to the matrix in (4.17). The extreme cases are Hqt(p) 



(0) and 



H 



Lp ] = (°)> while #imP is empty. 

r(*) 



We already know that Hy* p is strictly positive definite. As claimed earlier, its spectrum is 
R + and purely continuous, with a non-normalizable zero mode given by (4.18). Being finite- 
dimensional, the other two blocks jointly yield again just r non-negative and non-isolated proper 
eigenvalues for fixed fc>0. Interestingly, two of these eigenvalues are now zero and give us one 
even and one odd complex normalizable zero mode. With (3k, 7k £ C the latter take the following 
form: 18 



18 The even zero mode is of course absent for k — r. 
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H 



(k) 



imP ■ 



(j*e 



H, 



(*) 



Gr(P) • 



Ik 



0k 



■ • • y/r-k-2y/r-k-l\ 
v / FkL v / 2v / 3 • • ■ y/r-k-2y/r-k-l 
Vk+l^/k+2V3-- ■ y/ r -k-2y/r-k-l 



^k+ly/k+2^/k+3 ■ ■ ■ ^fr~^2^Jr-k-\ 
\Vk+lVk+2~Vk+3 ■ ■ ■ ^i^2^F^l 



( y/r+ly/r+2y/r+3 ■ ■ ■ Vr+k-2y/r+k-l \ 
^r-k+l^r+2^r+3 ■ ■ ■ ^r+k-2^r+k-l 
Vr-k+Wr-k+2^r+3 ■ ■ ■ Vr+k-2Vr+k-l 



y/r-k+Wr-k+2^r-k+3 ■ ■ ■ y/r-2^r+k-l 



\Vr-k+Wr-k+2y/r-k+3 ■ ■ ■ y/r^y/r-i j 



(4.24) 



(4.25) 



Altogether, there are 2r— 2 real zero modes in H im p and 2r real zero modes in H Gr ^ P y Identifying 
H(k)j = A-j,k+j = —A-k+j,j, the latter precisely agree with the BPS moduli found in (3.26). The 
former zero modes correspond to moduli outside the Grassmannian, thus generating nearby non- 
diagonal non-BPS solutions to the equation of motion. 



4.3.3 Diagonal perturbations 

We come to the diagonal (or radial) perturbations 0( O ) £ £0 ( see (4-7)). The transformation of 

00 

±EW[<S> r ,<j> {0) ] = Y,(m+l)(<P m+ i-<p m ) 2 -2r(4> 2 r + 4> 2 r _ 1 ) (4.26) 

m=0 

to a sum of squares is much easier than that of (4.11), but the result shows one minus in the 
signature: 

^£ (2) [*r,0(O)] = -r{^r-l+4>rf + £ (<f> j+1 - fa) 2 . (4.27) 

3>0,&r-l 

However, we can use this expression to conclude that the second variation form E^ cannot have 
more than one negative mode on So- 

Lemma. Let V C So be a real vector subspace such that [$> , <fi] < for all nonzero <f> G V. 
Then V is at most one- dimensional. In other words, one cannot find linearly independent vectors 
(f), Y> G £q such that E^[&, .] takes negative values on each non-zero linear combination of <fi and 

Proof. A necessary and sufficient condition for a pair of linearly independent vectors <f>, ip £ So 
to span a negative subspace for E^ 2 \.) := E^[^, .] is that 

E i2) (cp) < , E (2) (%b) < and B(cb,%b) 2 - E {2) {cp) E {2) (%b) < , (4.28) 

where B(<p, ip) := ^{E^ (<p+ip) — E^ (<f>) — E^ (Y>)} is the corresponding bilinear form. If 0, -0 G So 
satisfy these inequalities, then the same holds for their truncations (that is, vectors obtained by 
replacing all the coordinates (f>j, tpj with j > n by zero) if n is large enough. Thus, we may assume 
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the existence of n such that <f>j = ipj = for j > n. But the restriction of E^ to the space IR™ of 
all vectors <fi = (</>o, (pi, ... , 4> n -i, 0, 0, ... ) has signature (1, n— 1) according to (4.27) and, therefore, 



this space contains no two-dimensional negative subspaces. 

For an alternative visualization, we write the second-order fluctuation form as 



m,e=o 



(0) 

ml 



where the matrix entries can be extracted from (3.6) or from (4.27) as 

h H m£ = ( 2m + l ~ 2r<5 mjf ._i - 2r5 m)T ) 5 m/ - m 5 m/+1 - (m+1) S m+1: 
In matrix form this reads as (m,£ £ No) 



i ■ 



□ 



(4.29) 



(4.30) 



( 1 

-1 



V 



-1 

3 

-2 



-2 
5 
-3 



-3 
7 



- 2r(S. 



m,r—l' J £,r—l 



+ s n 



2r-3 
-r+1 



-r+1 
-1 

—r 



—r 

+1 -r-1 
-r-1 2r+3 



V 



(4.31) 
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In distinction to the earlier cases, this infinite Jacobian (i.e. tri-diagonal) matrix does not split. We 
shall compute its spectrum in a moment. Prior to this, some observations can already be made. 

First of all, the unique zero mode of is found from (4.18) for k=0 and given by (70 G R) 

4>zero = l70$r (</>m)zcro = 70 (~1, • • • , -1, +1, +1, ■ ■ ■) , (4.32) 



r times 

which clearly generates the global phase rotation symmetry $ — ► e 170 ^ already noted in (2.11). In 
contrast to the zero modes in £k>o depicted in (4.18), (4.24) and (4.25), zero is not Hilbert-Schmidt 
but bounded. Although it does not belong to a proper zero eigenvalue of the Hessian, 4> ze ro, being 
proportional to <3? r , still meets all our requirements for an admissible perturbation. 

In order to identify unstable modes we turn on all diagonal perturbations in a particular manner 
suggested by (3.30), 

) — 1 oo 

$({a m }) = 1 - ^2(e iam +l)\m)(m\ + ^ (e icVm -l) \m)(m\ . (4.33) 

m=0 m=r 



The energy of this configuration is easily calculated to be 



±E({a m }) = sin 



2 qp — a; 



_l_ . . . _|_ (r—1) sin 



2 Or 



+ r cos 



2 a r -_i_ar_ + (r+1) gin 2 a r -a r+1 



■ + ••• 

(4.34) 
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with a single cos 2 appearing only in the indicated place. Expanding up to second order in the real 
parameters a m and putting a m = for m > r, we find that 



ff(°> < 



(a -a 1 ) 2 + 2(ai-a 2 ) 2 + --- + (r-l)(a r -2-a r - 1 ) 2 -ra^ 1 < 0, (4.35) 



which determines a convex cone in the r-dimensional restricted fluctuation space. The most strongly 
negative mode is given by (a € R) 



= — iaP r 



to J neg 



-a (1^^,0,0,...) . 

r times 



(4.36) 



Comparison with (3.34) reveals that following this mode one arrives at $ = 1. More generally, 

(4> m ) = -a(0 1 _^ I 0, 1^^1,0,0,...) (4.37) 

f times (r—f) times 

perturbs in the direction of $ = 1 — 2Pf. None of these modes are eigenvectors of the matrix H^°\ 

It is instructive to confirm these assertions numerically. To this end, we truncate the matrix 
(H^j) at some cut-off value m max = £ max and diagonalize the resulting finite-dimensional matrix. 
Truncating at values m max < r only allows for positive eigenvalues, whereas for values m max > r 
we obtain exactly one negative eigenvalue. Its numerical value depends roughly linearly on r and 
converges very quickly with increasing cut-off as can be seen from the following graphs. 



i i I I I I I I I I I I i I I I 



i i l i i i i i l 



; i j I i i i i i i i j i i i i i i j i 



iiiiiiii 



^max 



Figure 1: Possible eigenvalues A versus the cut-off parameter m max for r = 5 and r = 20 



4.3.4 Spectrum of the Hessian 

The spectrum of the restriction = — ($ r A$ r )^ of the Hessian H to can actually be 
calculated by passing to a basis of Laguerre polynomials. The following considerations are based 
on spectral theory and require all vectors to be in fa, i.e. all operators to be Hilbert-Schmidt. 

Let us first consider £q. We use the following property of the noncommutative Laplacian — Ao: 
sending each basis vector \m) of TC to the m-th Laguerre polynomial L m (x), we get a unitary map 
U :H — > L 2 (R + , e~ x dx) such that UAqU^ 1 is just the operator of multiplication by x: 

xL m (x) = -mL m -i(x) + (2m+l)L m (x) - {m+l)L m+1 (x) , (4.38) 
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which is to be compared with (4.31). Hence, in the basis of the Laguerre polynomials, the eigenvalue 
equation H^\<p) = X\(f>) is rewritten in terms of / := U\(f>) as 



x/(x)-2r«/,L r _i)L r _i(x) + (/,L r )L r (x)) = Xf(x). 
Clearly, a function / G L 2 (R + , e~ x dx) satisfies (4.39) if and only if it is given by 

c L r _i(x) + c\L r {x) 



m = 



x — X 



(4.39) 



(4.40) 



and the constant coefficients cq, c\ G R satisfy the linear system 

(l r -l,r-l(X) - C + / r -l,r(A) Ci = 
i^r-^A) C + (ir,r(A) - Ci = 

which is obtained simply by inserting (4.40) into (4.39). Here we use the notation 



(4.41) 



x — A 



Lfc(x) Li(x) for k,l>0. 



(4.42) 



The integrals (4.42) are simple (though not elementary) special functions of A. Indeed, ioo(A) 
is a version of the integral logarithm: ioo(A) = — e~ A li(e A ) for all A < 0. On the other hand, using 
the recursion relations 

(k+l)L k+1 (x) = (2k+l-x) L k {x) - fcL fc _i(x) , (4.43) 
one can show by induction over k and / that all functions Ik,i(X) are expressed in terms of Ioo(\): 

hi(X) = A kl (X)I 00 (X) + B kl {\) , (4.44) 
where A^i and B k i are polynomials in A of degree at most k+l. Hence, the determinant 



F r (X) :-- 



Ir- 



l,r-l 



(A) 



2r 



r,r— 1 



(A) 



Ir-l,r(X) 
I r ,r(X) — 2^ 



(4.45) 



of the linear system (4.41) is a known special function of A, whose zeros A r on the negative semi- 
axis are precisely the negative eigenvalues of the Hessian operator H^[& r ] for the diagonal BPS 
background of rank r. 

One can now prove the existence of negative eigenvalues. By verifying that the real numbers 
F r (— oo ) := lim^_ > _ 00 F r (X) and F r (0) := lim^o- F r (X) have different signs, we can conclude that 
F r (X) possesses at least one zero on the negative semiaxis. For example, take r = 1. In this case 
we find 



Joi(A) = /io(A) = (l-A)Ioo(A) - 1 

Jn (A) = (1-A) 2 / 00 (A) + A-1 



2Fi(A) = -A-i-A 2 / 00 (A) . (4.46) 



Passing to the limit under the integral sign, we see that Fi(— oo) = ^ and ^1(0) = —\. This proves 
the existence of a negative eigenvalue of 
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Can H(°)[$> r ] also have non-negative eigenvalues? To disprove this possibility, we consider the 
eigenvalue equation (4.39) and show that it admits only the trivial solution f{x) = if A > 0. 
Indeed, for non-negative A the solution (4.40) is not square-integrable unless 

coL r ._i(A) + c\L r {X) = =>■ Co = cL r (A) and c\ = — cL r _i(A) (4.47) 

for some constant c. Therefore, the solution (4.40) is completely determined as 

/(*) = ^WLr-^-L , { X)L r{x) = c g Lfe(A)Lfe(x) (44g) 

fc=0 

via the Christoffel-Darboux formula. It follows that / is orthogonal to L r . Inserting (f,L r ) = 
into (4.39), we learn that the polynomial (x—X)f(x) is proportional to L r _i(x). But then the 
first equality in (4.48) demands that L r _i(A) = 0, constraining A. Feeding this into the second 
expression for / in (4.48) we see that the sum runs to r— 2 only, implying that / is orthogonal to 
L r _i as well. This finally simplifies the eigenvalue equation (4.39) to the one for Ao, namely 

xf(x) = Xf(x) , whence f(x) = , (4.49) 

as required. Hence, the non-negative part of the spectrum of is purely continuous. 

A similar analysis can be applied to the Hessian on £k>o- In this case, the appropriate polyno- 
mials are the normalized (generalized) Laguerre polynomials L^, which form an orthonormal basis 
for L 2 (R_|_, x k e~ x dx). One rediscovers the r proper eigenvalues as zeroes of the characteristic poly- 

(k) (k) 

nomial for H^ p © H^tpy accompanied by a continuous spectrum covering the positive semiaxis 

for H^2 P , as claimed earlier. 

If we do not care for resolving non- isolated eigenvalues, we may also infer the spectrum of 
the Hessian H from the classical Weyl theorem (see [30], Theorem IV. 5. 35 and [31], § XIII. 4, 
Example 3). 

Assertion. Let $ = <I> r be the diagonal BPS solution of rank r > 0. Then 

spectrum {H[$ r ]) = {A r } U [0, +oo) , (4.50) 

where A r is a negative eigenvalue of multiplicity 1 and [0, +oo) is the essential spectrum. 19 Fur- 
thermore, we have — 2r < A r < 0. 

Proof. Weyl's theorem asserts that the essential spectrum is preserved under compact pertur- 
bations. The essential spectrum of A (on all subspaces £&) is known to be equal to [0, +oo). (This 
follows, for example, from the explicit diagonalization of A& in terms of the Laguerre polynomials, 
as indicated above.) Since the perturbation ($ r A$ r ) is finite-dimensional, Weyl's theorem applies 
to show that the essential spectrum of H is also equal to [0, +oo). Although we have explicitly 
shown that is positive semi-definite for fc>0, the essential spectrum R + cannot exhaust the 
whole spectrum because the operator is not positive semi-definite on £q. (Indeed, (4.27) 

shows that the quadratic form is negative on many finite vectors.) Hence, must have 
negative eigenvalues of finite multiplicity. The above lemma implies that there can be only one 
such eigenvalue, and its multiplicity must be equal to 1. This proves (4.50) and the inequality 
A r < 0. The estimate — 2r < A r follows because the norm of the operator (<l? r A<3? r ) is equal to 2r, 
so H[<& r ] + 2rl^ > A is non-negative definite and hence cannot have negative eigenvalues. □ 

19 We recall that the essential spectrum of a self-adjoint operator is the whole spectrum minus all isolated eigenvalues 
of finite multiplicity. The multiplicity of an eigenvalue is the dimension of the corresponding eigenspace. 
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4.4 Results for non-diagonal U(l) BPS backgrounds 

Even though any non-diagonal BPS background $ can be reached from a diagonal one by a unitary 
transformation U which does not change the value of the energy, the fluctuation problem will get 
modified under such a transformation because 

Uatf = f{a,a)) and U a)U^ = f\a,a)) (4.51) 

will in general not have an action as simple as a and a) in the original oscillator basis. The exceptions 
are the symmetry transformations, of which only the rigid translations D(a) and rotations R("&) 
keep $ within the Grassmannian. To be specific, we recall that 

D{a) ] aD{a) = a + a =► E[D{a)$D{otf\ = E [$] , 

(4.52) 

R^ aR(#) = e^a => E [R(#) $ i?(tf) f ] = E [$] . 

The invariance of the spectrum of H under translations (or rotations) can also be seen directly: 
H4> n = e n 4> n and A(Dfrf) = D (A/) D* => 

(4.53) 

H'4>' n = [A-(D$D^A(D<S>D^]D<p n D^ = D ([A - ^>A^]<p n )D^ = e n <f>' n . 

More general unitary transformations will not simply commute with the action of A and thus 
change the spectrum of H. 

In the rank-one BPS situation, any solution is a translation of <E>i, and hence our previous 
discussion of fluctuations around diagonal U(l) BPS backgrounds completely covers that case. 
This is no longer true for higher-rank BPS backgrounds. For a complete stability analysis of 
abelian r-solitons it is therefore necessary to investigate separately the spectrum of H around each 
soliton configuration (2.56) based on r coherent states, whose center of mass may be chosen to be 
the origin. Since the background holomorphically depends on r parameters a±, . . . , a r and passes 
to <lv when all parameters go to zero, one may hope to show that the spectrum of H does not 
change qualitatively when $ varies inside the rank-r moduli space. 

5 Perturbations of U(2) backgrounds 

In this section we investigate the behavior of a simple U(2) BPS solution under perturbations and 
conclude that the fluctuation analysis reduces to the U(l) case discussed in the previous section. 

5.1 Results for diagonal U(2) BPS backgrounds 

We consider a nonabelian BPS projector of the diagonal type considered in (2.62): 

P = In © Pq where 
Clearly this describes a BPS solution 

$ ~ 1 $(2,1) $(2,2) j 



k=0 



-In 
1 H - 2P Q 



(5.1) 



(5.2) 
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of energy E = 8irQ. Inserting this into the expression (3.4) for the quadratic energy correction E^ 
one obtains 

£( 2 ) = 2 7 rTrU t A0 - $ (*A$t) <f>] 

(5.3) 

= 27rTr{</> t U?)A</> - (g ?) 2Q(|Q-1)(Q-1| + |Q)(Q|) </>} • 
Since the action of the Hessian is obviously diagonal, we find that 

which reduces the fluctuation analysis to a collection of abelian cases. For the case at hand, the 
Hessian in the (1, .) sectors is given by the Laplacian and thus non-negative, while in each (2, .) 
sector it is identical to the Hessian for the abelian rank-Q diagonal BPS case. Hence, the relevant 
results of the previous section carry over completely. 

This mechanism extends to any diagonal U(2) background (not necessarily BPS), mapping the 
U(2) fluctuation spectrum to a collection of fluctuation spectra around corresponding diagonal 
U(l) configurations. The fluctuation problem around non-diagonal backgrounds, in contrast, is not 
easily reduced to an abelian one, except when the background is related to a diagonal one by a 
rigid symmetry, as defined in (2.11), (2.13) and (2.16). Only a small part of the moduli, however, 
is generated by rigid symmetries, as our example of U(fi) in (2.68) demonstrates even for Q=l. 
Finally, it is straightforward to extend these considerations to the general U(n) case. 

6 Conclusions 

After developing a unified description of abelian and nonabelian multi-solitons in noncommutative 
Euclidean two-dimensional sigma models with U(n) or Grassmannian target space, we have ana- 
lyzed the issue of their stability. Thanks to the BPS bound, multi-solitons in Grassmannian sigma 
models are always stable. As in the commutative case, their imbedding into a unitary sigma model 
renders them unstable however, as there always exists one negative eigenvalue of the Hessian which 
triggers a decay to the vacuum configuration. 

Our results are concrete and complete for abelian and nonabelian Q-soliton configurations which 
are diagonal in the oscillator basis (or related to such by global symmetry). For this case, we proved 
that the spectrum of the Hessian consists of the essential spectrum [0, oo) and an eigenvalue Xq 
of multiplicity one with — 2Q < Xq < 0. This assertion was confirmed numerically, and the value 
of Xq was given as a zero of a particular function composed of monomials in A and the special 
function e _A li(e A ). 

Furthermore, the complete set of zero modes of the Hessian was identified. Each abelian diagonal 
Q-soliton background is characterized by a diagonal projector P of rank Q, whose image and kernel 
trigger a decomposition of the space of fluctuations into three invariant subspaces, namely u(imP), 
-u(kerP) and dGr(P). In addition, every side diagonal together with its transpose is seperately 
invariant under the action of the Hessian. This leads to a particular distribution of the admissible 
zero modes of the Hessian, displayed here for the example of Q=4: 
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where the double line denotes the single negative eigenvector, each solid diagonal segment represents 
a real normalizable zero mode, the dashed line depicts an admissible non-normalizable zero mode, 
and empty areas do not contain admissible zero modes. In addition, each side diagonal features a 
non-admissible zero mode at the edge of the continuous part [0, oo) of the spectrum. We plot the 
complete spectrum of the Hessian at Q=A (cut off at size m max =30) for each invariant subspace 

S^ 1 ^ (boxes), £™ p (stars), £^ crP (crosses) and £q (circles), up to k=6: 



15 - 
12.5 - 
10 - 
7.5 - 
5 - 
2.5 - 



Figure 2: Eigenvalues A of the cut-off Hessian (size 30) for Q=4 in subspaces 



Since most soliton moduli are not associated with global symmetries, our explicit results for 
diagonal backgrounds do not obviously extend to generic (non-diagonal) backgrounds. We have 
not been able to compute the fluctuation spectrum across the entire soliton moduli space. The 
only exception is the abelian single-soliton solution which always is a translation of the diagonal 
configuration and therefore fully covered by our analysis. Already for the case of two U(l) solitons, 
the unitary transformation which changes their distance in the noncommutative plane is only 
partially known. 

This leads us to a number of unsolved problems. The most pressing one seems to be the extension 
of our fluctuation analysis to non-diagonal backgrounds. Next, following the even zero modes one 
can now find new non-BPS solutions. Also, it is worthwhile to investigate the commutative limit of 
the Hessian and its spectrum. Another interesting aspect is the existence of infinite-rank abelian 
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projectors, i.e. via an infinite array of coherent states, associated with BPS solutions of infinite 
energy. Furthermore, some technical questions concerning the admissible set of fields and their 
fluctuations have remained. Finally, it would be rewarding to extend the geometrical understanding 
of sigma models to the noncommutative realm. 
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